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Abstract 



An algebraic approach for factorizing nonlinear partial differential equations (PDEs) 
and systems of PDEs is provided. In the particular case of second order linear and non- 
linear PDEs and systems of PDEs, necessary and sufficient conditions of factorization 
are given. 
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1 Introduction 

The search for exact solutions of differential equations is very challenging in mathematics, 
but their usefulness in the proper understanding of qualitative features of phenomena and 
processes in various areas of natural science merits to get down to such an investigation. 
Indeed, exact solutions can be used to verify the consistency and estimate errors of various 
numerical, asymptotic and approximate analytical methods. Unfortunately, there does not 



*E-mail address: norbert.hounkonnou@cipma.uac.bj 
'E-mail address: sielenou_alain@yahoo.fr 



always exist a method adapted for the resolution of any type of differential equations. Very 
often, one tries to reduce the equation in order to make easier its resolution. But this reduc- 
tion requires the knowledge of suitable transformations or changes of variables. The latters 
usually give rise to another problem the issue of which is not always favourable. 

A simple approach for the reduction of a differential equation consists in seeking a 
factorization, if there exists, of the differential operator associated with it. Note that for 
the particular case of second order linear ordinary differential equations of Schrodinger or 
Sturm-Liouville type, the factorization of the associated differential operators also allows 
to obtain partially or completely their spectrum, under certain assumptions of integrability 
|[T]|2l[3l|4l. In recent years, there has been much interest devoted to the problem of factor- 
ization of differential equations, especially based on linear ordinary ['5^, '6', '8^ and nonlinear 
differential operators tZ,^,^. Although effective, the used methods are rather restrictive in 
their applications. 

Recently, a purely algebraic method of factorization of the second order linear ordinary 
differential equations has been presented by the authors in lilOl [TTl IT2l [131. The same 
procedure of factorization has been exploited in |[T4l and extended to second order nonlinear 
ordinary differential equations (NLODES) and systems of NLODEs. This work generalizes 
previous works by applying the above mentioned algebraic method of factorization to linear 
and nonlinear systems of partial differential equations (PDEs). Necessary and sufficient 
conditions of factorization are derived in the case of second order equations. 

First of all, some useful notations are required. Consider X, an n-dimensional in- 
dependent variable space, and U, an m-dimensional dependent variable space. Let x = 
(x^ • • • ,x") e X and M = {u\--- ,u'") E U. We define the space [/(') , s G N as: 



^W:= „W:„W = (^ (^J I (1.1) 




where uLs is the 



Pk=n'' (1.2) 



of all k-th order partial derivatives of u-' . The M/^n vector components are recursively ob- 
tained as follows: 

i) "(0) = uJ and m^^^ = (^m^i , m^j , • • • , m^, j ; 
ii) Assume that ujj^. is known. Then, 

- Form the tuples jH^ ^^ (Z) as follows: 

where M/^v [/] is the Z-th component of the vector M/^^ ; 

- Finally, form the vector 

^ -'"^ i)(l)'«(/t+i)(2)r--,"Ln(w) 



Un , ,\= U 



{k+i)~ y'{k+\)y'-)->'^[k+\)y^>-> ■>'^{k+\ 
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An element u^'^\ in the space U^'^\ is the 

qs = m{l+pi+p2-\ h As)-tuple (1.3) 

defined by 

The coordinates in the space X x U^-'^^ are denoted by (jc, m^^)) . 

In the sequel, the qs-uple u^^' will be refen^ed to (I1.4I ). whereas the integers pk and q^ are 
defined by (11.21 ) and (11.31 ). respectively. Define differential operators Dt^h whose action on 
a regular function u is 

DkjrU = u^k)[h] (1.5) 

These operators Djt,/, satisfy the following properties: 

(i) Do,i u = u (identity), 

(ii) Di ftD^/ft/M = D|./_|_i „(/,/_i)^/, M (composition rule), 

(iii) Djt/jM = Di_/, Di_i 1 M, k>l (decomposition rule). 

Remark 1.1. Operators D^/, allow the simplification of the writing of certain differential 
operators. For example, the operator 



can be shortly expressed as 

.V Pk 



2 Linear differential operators 

In this section, we develop an algebraic method of factorization applicable to linear differ- 
ential operators (LDOs) and to systems of LDOs. 



2.1 Factorizations of linear differential equations 

The general setting of the factorization problem for LDOs is developed. Necessary and 
sufficient conditions are derived for the factorization of second order linear ordinary and 
partial differential operators with two independent variables. 
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2.1.1 General setting 

Let 5 > 2 be a positive integer and A be an open subset of M." . Let 

s Pk 

2'(^) = ££guWDu (2-1) 

/t=0/!=l 

a linear differential operator of order s, where gk^h £ C (A, M) . The operator T [s) acts on a 
function m G (r*(A,M) as follows 

i Pk 
'P{s)u=Y^Y^ gk,h ix)Dkji u. (2.2) 

k=Qh=l 

The method of factorization consists in seeking a decomposition of the differential operator 

(12.11) in the following form 

/ 
'P{s)=Y\ch{si) (2.3) 

i=\ 

with ^-^1 Si = s and 



Si Pk 

II 

k=Oh=l 



Qj{si) = 11 ba,h{x)V)k.h, (2.4) 



where Z7i,jt,/j G (:(A,M) and bi^k.h G c'^'j='''{A,R), i = 2,3,--- ,1. 

Proposition 2.1. Let P{s) be an operator which can be decomposed into the form d2. JD . If 
the function uq satisfies 

Ql{si)uQ = 0, (2.5) 

and Ml , . . . , M;_i are solutions of the system 

I 
n Qjc{^k)uj = Vj, j = l,2,...,l-l, (2.6) 

k=i-j+i 

where Vy, 7 = 1 , 2, . . . , Z — 1 , are solutions of 

l-j 

Y[ai{si)vj = 0, (2.7) 

f/jen Mo, Ml , . . . , M/_i are I particular solutions of the equation 'P[s)u = 0. 
Proof. Let mq and mj, j = 1, 2, . . . , Z — 1 be solutions of (12.51 ) and (12.61) . respectively. Then 

//-I \ 



^(^)"o= nQ*' (■*'■) Q/('5/)"o = 0, 



and for j = 1, 2,..., Z — 1, 

'P{s)uj = \Yl(h{si)\{ n Clk{sk)\uj 

l-j 

(=1 
where the use of (12.61) and (12.71) has been made. D 
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Expanding (12.31 ) leads to the relations between unknown functions Z?;,^,/, of the differen- 
tial operators Qi{si) and the known functions g/t./j of the original differential operator 'P{s). 
Without loss of generality and as matter of clarity, this study will be concentrated to second 
order equations, the generalization being straightforward. 



2.1.2 Necessary and sufficient conditions for the factorization of second order linear 
ODEs 

Let A and Aq be two open subsets of M such that Aq C A. Consider the second order linear 
ordinary differential operator 

2 Pk 
2' (2) = Y.l.Sk,h{x)^Kh 

k=Oh=\ 
= ^0,1 WDo,l +gLl(x)Di4 +^2,1 WD2,1, 

where gk.h G C (A, M) and x = xK Write (P (2) in the form 



(2.8) 



2'(2) = Q,i(l)-Q2(l^ 
1 Pk 



k=Oh=l 



1 Pk 
Y, X^21,/!WD<:,/, 
k=Oh=l 



(2.9) 



(2.10) 



= [^lAl W^O,! +^1,1,1 WDl.l] [^2,0,1 WDo,l + ^2,1,1 WDli] : 

where bi^t,h S (:(A,M) and bj^kji G (:^(A,M). Let u £ C^{Ao,M.). Then we have 

f{2)u = gOA{x)u+gi,j{x)Ux+g2A{x)u2x 

and after expansion 

2' (2) M = [bifl^i (x)Do,i +biA,i WDi,i] [^2,0,1 WDo,i + ^2,1,1 WDi,i] " 

= ^1,1,1^2,1,1 U2x + [^1,0,1^2,1,1 + ^1,1,1^2,0,1 +^i,i,iDi,i (^2,1,1 )] Ux 

+ [^1,0,1^2,0,1 +^1, i,iDi, 1(^2,0,1)] "• (2.11) 

Identifying dTTOl ) with drTTT ) yields 

Proposition 2.2. A necessary and sufficient condition for the differential operator (B (2) 
defined by f |2.8D be decomposed into the form f |2.9D is: 



glA = ^1,1,1^2,1,1, 

glA = ^1,0,1^2,1,1+^1,1,1^2,0,1+^1,1,1^1,1(^2,1,1), 

^0,1 = ^1,0,1^2,0,1+^1,1,1^1,1(^2,0,1)- 



(2.12) 
(2.13) 
(2.14) 



Propose an approach to solve system (12. 12I )- (I2. 141 ). Assume that g2.\ does not vanish on 
A. Thus, it is always possible to find two nonzero functions on A, namely ^1.1,1 and ^2,1,1, 
which satisfy (12.121 ). Substituting X = b\^Q\ and Y = ^2,0,1 in (12.131) gives 



X 



1 



b2. 



1,1 



[gi,i -^1,1, iDi, 1(^2,1,1) -^1,1,15^] • 



(2.15) 
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The substitution of (12.151) into (12.141 ) implies that the decomposition ( 12.91 ) is strongly related 
to the existence of a solution to the following Riccati equation in Y 



DM(n 



g2A 



82,i 



gQ,l 



'1,1,1 



0. 



(2.16) 



2.1.3 Necessary and sufficient conditions for the factorization of second order linear 
PDEs with two independent variables 

Let A and Aq be two open subsets of M? such that Aq C A. Consider the second order linear 
partial differential operator 

2 Pk 

2' (2) = Y.LskA^)^Kh 

k=Qh=l 
= ^0,1 WDo,l +^Ll(x)Di4 +^L2WDi^2 
+ g2,lix)'D2A+g22{x)D22+g23{x)D23+g2A{x)D2,A, (2.17) 

where gk.h G (r(A,M) and x = (x' ,x^) . Write !P (2) in the form 



2'(2) = Qj(l)-Q2(l) 



1 Pk 
k=Qh=l 



1 Pk 

_k=Oh=l 

= [^1,0,1 WDo,i +^l,l,lWDl,l +Zji,i.2WDi,2] 
X [^2,0,lWDo,l+Z72,l,l(x)Di,i+Z;2,l,2WDi,2], (2.18) 

where b^^kj, £ c{A,R) and b2,k,h e C^{A,R). Let u £ C^{Ao,R). Then we have 

(P{2)u = go,l" + gl,l"xl +gl,2'^x^+g2,lU2x> +(52,2+52,3) U^ijc2+g2AU2x^i2.l9) 

and after expansion 

'P{2)u = [^i,o,iWDo,i+^i,i,iWDi,i+^i,i,2WDi,2] 
X [^2,0,1 WDo,i + ^2,1,1 WDi,i + ^2,1,2 WD!,2] " 

= [^1,0,1^2,0,1 +^1, l,lDl, 1(^2,0,1) +^l,l,2Di,2(^2,0,l)] U 

+ [^1,0,1^2,1, 1+^1,1, 1^2,0,1 +^l,l,lDij(Z72,l,l)+^l.l,2Dl,2(^2,l,l)] ^x' 

+ [^1,0,1^2,1,2 +^l,l,2%),l+^l,l,lDij(ft2,l,2) +^1,1,201,2(^2,1,2)] "a:2 

+ ^1,1,1^2,1,1 "2;f' + [^1,1,2^2,1,1+^1,1,1^2,1,2] "xl.v2 +^1, 1,2^2,1,2 «2x2- (2-20) 

Identifying ( 12.191 ) with ( 12.201 ) leads to the following 

Proposition 2.3. A necessary and sufficient condition for the differential operator V (2) 
defined by ([2. 1 71 ) be decomposed into the form d2.i8D is: 

52,1 = ^1,1,1^2,1,1, (2.21) 

52,2+52,3 = ^1,1,2^2,1,1+^1,1,1^2,1,2, (2.22) 

52,4 = ^1,1,2^2,1,2, (2.23) 

51.1 = ^1,0,1^2,1,1 +^1,1,1^2,0,1 +i^(^2,i,i), (2.24) 

51.2 = ^1,0,1^2,1,2+^1,1,2^2,0,1+^^(^2,1,2), (2.25) 
50,1 = ^1,0,1^2,0,1 +-^(^2,0,2), (2.26) 
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where L = ^i,i,iZ)i i +Z7ij 2^1 2. 

Propose an approach to solve system (12.21 1 )- (I2.26I ). Assume that at least one of the 
functions g2.\ and ^2,4 does not vanish on A, says ^2,1- It is always possible to find two 
nonzero functions on A, namely ^1.1,1 and ^2,1,1 which satisfy (I2.21I ). Substituting X\ = 
bi,\^2 and X2 = ^2.1,2 into (12.221 ) yields 

^1 = 7 (52,2+^2,3 -^1,1,1^2). (2.27) 

»2,1,1 

The substitution of ( |2.27b into ( 12.231 ) shows that ^^2,1,2 is a solution of the second degree 
algebraic equation 

^W ^2 _ g2,2+g2,3 ^^ + 1^=0. (2.28) 

g2.l g2A ^1,1,1 

The discriminant of equation (12.281 ) is 

A = (^2,2 + 52,3) -4g2,i§2,4 = (^1,1,2^2,1,1 -^1,1,1^2,1,2) >0. (2.29) 

If A > 0, then the substitution of F = ^1,0,1 and Z = ^2,0,1 into ( 12.241 ) and ( 12.251) implies 
that the decomposition (12.181 ) is possible if the unique solution to the following algebraic 
system in Y and Z 

gl,l-L{b2,l,l) = ^2,1,1^^ + ^1,1,12 

§1,2-^ (^2,1,2) = ^2,1,2^^ + ^1,1,22 (2.30) 

satisfies ( 12.261 ). Indeed, the determinant of the system ( 12.301 ) is 

^1,1,2^2,1,1 - ^1,1,1^2,1,2 = ±Va / 0. 

If A = 0, then the substitution ofY = ^1,0,1 and Z = ^2,0,1 into (12.241 ) yields 

Y = ^ ki,i - L (^2.1,1) - ^i,i,iZ] . (2.31) 

»2,1,1 

Then, the substitution of (12.311 ) into ( 12.261 ) implies that the decomposition (12.181 ) is strongly 
related to the existence of a solution to the following first order quasi-linear partial differ- 
ential equation in Z 

X (Z) - h^z' + Sui-^ihu) ^ -^01=0 (2.32) 

^2,1,1 ^2,1,1 

which satisfies (12.251 ). 



2.2 Factorizations of systems of linear differential equations 

The previous analysis is now made for systems of linear differential equations. 
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2.2.1 General considerations 

Let A be an open subset of M". Examine now the factorization process for systems of s-th 
order, {s>2), linear differential equations with n independent variables x= (x^ ,■■■ ,x"^ and 
m>2 dependent variables u = ' (m^, • • • ,m'") , u = u{x) whose associated matrix operator, 
fW (s) , is of the form 

5i/(.) = [%,(.p,,)]i<^,^<,„; (2.33) 

the lip^q i^p,q) ^e ■^p.^'th order linear differential operators 

^P-l Pk 

^,q (.^P,q) = L L fP'1-l<J^ (^)^'t.'' ' (2.34) 

k=()h=[ 

where fp,q,k,h G C (A, M) , Sp^g = ^ - 1 + 5^,^, 5^.^ = 1 and bp^q = Oif pj^q. 

Let A and Aq be two open subsets of M" such that Aq C A. The matrix operator !M (s) acts 

on a vector valued function u= ' [u^ ,■■■ ,u'"^ € C*(Ao,M'") as follows 



2^ ^,q \^p,q} 
q=\ 



— 1 

l<p<in 

The method of factorization consists in seeking a decomposition of the matrix M {s) under 

the following form 

/ 
M {s) =Y{9{i{si) (2.35) 

r=l 

where 

5V^-(^/) = [T,-,p,,(.,-,p,,)]i<^^^<^ (2.36) 

and 

^i.P'l Pk 
Ti.p,q {Si,p,q) =YjH <^i,p,qMA^)^k,h, (2.37) 

k=Qh=i 

with l!i=\Si = s, Si^p^q = Si - 1 + 5p,^, ax^p^q^k^h G (:(A,M) and ai^p^q^k,h G C^^=i''"'(A,R), 
/ = 2,3,-- • ,/. 

Proposition 2.4. Let M {s) be a matrix of differential operators defined by d2. JJI ) which 
can be decomposed into the form d2.J5D . If the function uq = ' (mq, • " " i^o ) satisfies 

^{s,)uQ = 0, (2.38) 

and Uj = ' iup-- ■ ,ujj , j = \,2,-- ■ ,1 — I are solutions of the system 

1 
n ^i^k)uj = Vj,j = \,2,...,l-\, (2.39) 

k=i-j+i 

where v^- = ' ( vj , • • • , v'" j , 7 = 1 , 2, . . . , Z — 1 , are solutions of 

i-J 

llo^{si)vj = 0, (2.40) 

1=1 

then uq, ui, .. . , ui-i are I particular solutions of the equation M {s)u= 0. 



Factorization of differential operators 



Proof. The proof is similar to that of the Proposition 12.11 D 

Expanding (12.351 ) leads to the relations between the unknown functions at^p^q^tji of 
!y{i{si) and the known functions fp^q,k,h of ^ [s). 

As matter of clarity, in the sequel we explicitly derive necessary and sufficient condi- 
tions for the factorization of systems of second order linear ordinary and partial differential 
operators with two independent variables. 

2.2.2 Necessary and sufficient conditions for the factorization of systems of second 
order linear ODEs 

Let A and Ao be two open subsets of R such that Aq C A. Consider the matrix operator 

^i^) = K,,U<p,q<,n^ (2-41) 

where 

2 Pk 

^,p = X, IL fp,g,k,hi^)^k,h = fp,p,Q,l^Q,l +/p,p,l,lDl,l + fp,p,2,l^2,\ (2.42) 

k=Oh=l 

and for p j^q 

1 Pk 
^,? = LI IL /p,?,M W^Wi = fp,qfl,i'DQ,i +/p,?,i,iDli (2.43) 

k=Qh=l 

with fp,cj,k.h £ CiA,R),x = x\ Write fW (2) in the form 

fW-(2) = 5V:i(l)-fA6(l), (2.44) 

where 

^•(l) = [T,-p,,]i<p^^<„ (2.45) 

with 

1 Pk 

'^',p,p ~ L L ^i:P:P,k,h{x)f^k,h = ai.,p.,pfi,\f^Q,\ +ai,p,p,\,\f^\,\ (2.46) 

k=Oh=Y 

and for p j^q 

"Ti^pxi = «!,p,?,o,iDo,i , (2.47) 

a\,p,q,k,h G (:(A,M) and a2,p,?,/t,/i G C^AjM). Let m = * {u\--- ,u'") G C^{Aq,R). Then we 
have 



where 



L^,?"'^ 
.«=i 



l<p<m 
^),pU'^ = fp,p,0,l u'^ +fp,p,i,\ K+fp,p.2,l "2x 



(2.48) 



and for p y^q 

'Kp,qU'' = fp,qfl,\ u'^+fp,qA,l M^ 



10 
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On the other hand, after expansion of (12.441) . we have 



fW(2)M 



'Kp,q 



\<p,q<m 






(2.49) 



\<p<m 



where 



^j.pU'' 



2^ ^l,p,/,0,lfl^2,/,p,0,l + ^l,p,p,l,lDl,l {a2,p,p,Q,l , 



u'^ + <^l.p.p.l,l<^2,p,p,l,l "2x 



+ [«l,p,p,0,l«2,p,p,l,l +fl'l,p,p,l,ia2,p,p,0,l +fl:i,p,p,l,lDl,l('32,p,p,l,l)]M^ 

and for p / ^ 

^.qU^ = [ai,p,p,l.ia2.p,q,0,l +ai^p^qfija2,q^q^l,l]uf 



+ 



}_i<^l,p,lfi,l<^2,l,q,0,l +«l,p,p,l,lDl,l(«2,p,?,0,l) 



/=1 



Identifying (12.481) with ( 12.491 ) yields 

Proposition 2.5. A necessary and sufficient condition for the differential operator fW (2) 
defined by Ii2.41i be decomposed into the form Ii2.44i is: 



fp,P,o,\ — 2j^i'P'''0,ifl:2,/,p,o,i +«i,p,p,i, 1-01,1(0^2,^,^,0,1)) 



(2.50) 



i=\ 



fp,p,l,l — '^l,p,p,0,ia2,p,p,l.l +ai,p,p,l.ia2.p.p,0,l +ai.p.p.l, 1^1,1 (a2,p,p.l,l), (2.51) 

fp,P,2,i = «i,p,p,i,i«2,p,p,i,i (2.52) 



and for p ^ q 



fp,qfi,l — 2^«l,p,/,0,lfl^2,/,^,0,l +ai.p.p.l,l^l,l(a2,/7,<?,0,l! 



Jp,q,l, 



1=1 

fl^l,p,p,l,l'32,p,fl,0,l +«l,p,o,0,l«2,fl,o,I,l- 



(2.53) 
(2.54) 



2.2.3 Necessary and sufficient conditions for the factorization of systems of second 
order linear PDEs witii two independent variables 

Let A and Aq be two open subsets of R^ such that Aq C A. Consider the matrix operator 



^(2) = [%,]i<p 



,q<m ' 



(2.55) 



where 



2 Pk 

k=Oh=l 
= /p,p,0,lDo,l +fp.p.l,l^l,l +fp,p,l,2^1,2 
+ /p,p,2,lD2,l +/p,p,2,2D2,2 +/p,p,2,3D2,3 + fp,p,2A^2,4 



(2.56) 
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11 



and fov p^q 



1 Pk 



^,q = 2_i 2_^ fp,q,k,h(^)^k,h = fp,q.Q,1^0,l + fp,q,l,l^l.l + fp,q,l,2^l.2 
k=Oh=l 

with fp^cj,k,h S (:(A,R), X = {x^,x^) . Write fW (2) in the form 

^(2)=fA£i(l)-5V2(l), 
where 



(2.57) 



^{i) = [Ti,p,q], 



<p,q<m 



(2.58) 
(2.59) 



with 



1 Pk 



"^'■p-p ~ H H ai,p,p,k,h{^)^k,h = «i,p,p,o,iDo,i +«(,p,p,i,iDi,i +«i>,p,i,2Di,2 (2.60) 



and for p^ q 



'^i,p,q —<^i,p,q,0,l^O,l, 



^i.p.q.kM G C(A,M) and a2^p^q^k,h £ C (A,M). Let u = ' [u 
have 



f ''"1 •••,«'" 



(2.61) 

) G C^(Ao,M).Then we 



^(2)" = fe<?]i<p,,<« 






(2.62) 



l<p<m 



where 



'Kp,pU^ — f p,p,0,luP +fp,p,l,l U^i+ fpp i^2t^^2 

+ fp,p.2,\ "2x1 + ifp,p,2.2 + fp,p,2,3 ) "^lj^2 + fp,p,2A "2;c2 

and for /? / ^ 

'Kp,qU'^ = fp,qfi,llfl + fp.q,l,lU^l + //;,5,1,2 "^2 ■ 

On the other hand, after expansion of (12.581) . we have 



M{2)u 



'Kp,q 



l<p.q<m 



£%^"'^ 



(2.63) 



\<p<m 



where 



%p,pU'' = <^l,p,p,l,l<^2,p,p,l,l "2jc1 +('^1.P>P>1,2'^2,p,p,1,1 + <^l,p,p,l,l<^2,p,p,l,2) U^ 

+ [<^\,p,p,Q,i<^2,p,p,\,i +a\,p,p,i,ia2,p,pfl,i +-^p(«2,p,p,i.i)]"^i 

+ [«l,p,p,0,l'^2,p,p,l,2+«l,p,p,l,2fl'2,p,p,0,l +-^p('^2,p,p,l,2)]"^2 



+ 



2^a:i,p,/,o,i'^2,/,p,o,i +^p(«2,p,p,0,i' 
1=1 



U^ + «l,p,p,l,2«2,p,p,l,2 "2a:2 
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and fov p^q 



'Kp,qU'' — [a\,p,p,\^a2,p,qfi^+a\,p,qfi^Cl2,q,q,\^]u^l 
+ V^l,p,p,l,2<^2,p,qfi,l + a\,p,qfi,l<^2,q,q,l,2]^x^ 



+ 



2_i ai,p,lfi,\a2,l,qfi,l + Lp{a2,p,qfl,\ ] 



1=1 



u\ 







Jp,p 


1.1 






fp,p 


1,2 






Jp,p 


2,1 




fp,P,22 


+ fp^P 


2,3 






Jp,p 


2,4 


andfof 


P^q 







where Lp = ai,p,p,i,iDij +ai,p,p,i,2Di,2- From the Identification of (12.621) with (12.631) re- 
sults 

Proposition 2.6. A necessary and sufficient condition for the differential operator fW" (2) 
defined by d2.55D be decomposed into the form ( 12.581 ) is: 

m 

fp.p.0,1 = 2_^aipifija2^i,pfl^i+Lp{a2,p,pfl,i), (2.64) 

1=1 

<^\.p.pX),\<^2,p,p,\,i +«i,p,p,i,i«2,p,p,o,i + J^p{a2,p,p,\.\), (2.65) 

fl:ip,p,0,lfl'2,p,p,l,2 +fl^l,p,p,l,2fl'2,p,p,0,l +^p(«2,p,p,l,2), (2.66) 

«i,p,p,i,ifl'2,p,p,i,i, (2.67) 

fl:i,p,p,i^2fl'2,p,p,i,i +fl:i,p,p,i,ifl:2,p,p,i,2, (2.68) 

«l,p,p,l,2«2,p,p,l,2 (2.69) 



fp,q,0,l — 2_iai^pJfi^ia2J,qfi,]_+Lp{a2^p^q,Q^l), (2.70) 

/=1 

fp,q,l,l = <^l,p,p,l,l<^2,p,q,0,l +<^l.p,q,0,l<^2,q,q,l,l, (2.71) 

fp,q,1.2 = a:i,p,p,l,2«2,p,^,0,l +«l,p,?,0,l«2,^,?,l,2- (2.72) 



3 Nonlinear differential operators 

In this section, we investigate the factorization of nonlinear differential operators. 

3.1 Factorizations of nonlinear differential equations 

We start with general considerations and then deduce the main results on conditions of 
factorization. 

3.1.1 General setting and results 

Let s > 2 be a positive integer, A be an open subset of M" and D. an open subset of M. Let 

s Pk 

^i') = LL8kAx,-)^k,h (3.1) 

k=Oh=l 
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be a nonlinear differential operator of order s, where g^j, G C (A x n,M). The operator (P (s) 
acts on a function u G C^{A,Q.) as follows 



.V Pk 

k=Oh=l 



(3.2) 



The method of factorization consists in seeking a decomposition of the differential operator 
(13.11 ) in the following form 



'P{s)=X{Qi{si) 



i=\ 



with Y!i^i Si = s and 



Si Pk 



Qdi^i) = YjH bi^k,h{x,-)^k,h, 



(3.3) 



(3.4) 



k=Qh=l 



where Z7i^;t,/! e C{Axa,M.) and bi^k,h G C^^=i'^(A x fl,M), / = 2,3,-- • ,/. 

Expanding (13.31 ) leads to the relations between unknown functions bikj^ of the differential 

operators Qj{si) and the known functions gk^h of the original differential operator 2' (s). 

3.1.2 Necessary and sufficient conditions for the factorization of second order non- 
linear ODEs 

Let Q.,A and Aq be three open subsets of R such that Aq C A. Consider the second order 
nonlinear ordinary differential operator 

2 Pk 

2' (2) = ££g^,/,(x,-)D^,/, 

k=Oh=l 



= go,i(-^,-)Do,i+.?Li(-^r)Di,i+g2,i(-^,-)D2,i, 
where gkfi e (T (A x H, M) and x = x\ Write f (2) in the form 



(3.5) 



2'(2) = Q,i(l)-Q2(l) 



1 Pk 

k=Oh=l 



1 Pk 

X J^l^2,k,hi^,-)'Dk,h 
k=Oh=l 



= [bifi^i (x, •)Do,i +^Li,i (x, ■)Dia] [b2fi,i {x, ODq,! + ^2,1,1 (x, ■)Gi,i] , (3.6) 
where Z?i,i,/, £ c{AxD.,R) andb2i,h G c\Ax D.,R). Let u G C^{Aq,D.). Then we have 

'P(2)u = gQ^i{x,u)u + gi^i{x,u)ux + g2,i{x,u)u2x (3.7) 

and after expansion 

2' (2) M = [^^1,0,1 {x, ^Doj + ^1,1,1 (x, ODli] [^2.0,1 (-^^ ODo,! + ^2,1,1 {x, ODi.i] M 

= [^1,0,1^2,1,1+^1,1,1^2,0,1 +^1,1,iDi,i(^2,1,i)+^1,1,iDi,2(^2,0,i)m]";c (3.8) 
+ ^l,l,lDi,2(^2,l,l)"x + [^1,0,1^2,0,1 +^1,1, lDl,l(^2,0,l)] " + ^1,1, 1^2,1,1 "2a:- 



Identifying (13.71 ) with (13.8b furnishes 
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Proposition 3.1. A necessary and sufficient condition for the differential operator 2" (2) 
defined by i\3.5\) be decomposed into the form 0.61) is: 

(3.9) 
(3.10) 
(3.11) 
(3.12) 



g2,l = ^1,1,1^2,1,1, 

^1,1 = ^1,0,1^2.1, 1+^1, 1,1^2,0,1 +^1,1, 1^1,1(^2,1,1) +^1,1, 1^1,2(^2,0,1)", 
^1,1,1^1,2(^2.1,1), 







§0,1 



^1,0,1^2.0,1+^1, 1,1^1, l(^2,0,l)- 



3.1.3 Necessary and sufficient conditions for the factorization of second order non- 
linear PDEs with two independent variables 

Let A and Aq be two open subsets of M^ such that Aq C A. Let Q. be an open subset of M. 
Consider the second order nonlinear partial differential operator 



Til) 



Pk 



= L L^Wi(-^,-)D/t,A 

k=Oh=\ 
= ^0,1 {x, ODci +§1,1 {x, ODli +g\2{x, 001,2 + g2,l {x, •)D2,1 

+ g2,2{x, 002,2 +§2,3 (JC, 002,3 +§2,4(-^, 002,4, 



(3.13) 



where gi,/, G c{h>iQ., 
2' (2) = 



andx : 



ix\x^) 



Write fP (2) in the form 



Q,i(l)-Q2(l) 
1 Pk 
Y^ Y^bi^k,h{x,-)Dk.h 

k=Oh=\ 



1 Pk 
Y^ X1^2,u(-X,0O<:,ft 
yt=0/7=l 



[Z?i,o,l(x,OOo,l +Z7ij.i(x,0Oi,i +^l,l.2(x,0Oi,2] 
[^2,0,1 {x, OOo,l + ^2,1,1 {x, OOl.l + ^2,1,2 (JC, 0Ol,2] 



(3.14) 



where bi^h £ (:(A x fl,R) and b2.k.h G C^^ x ^,IK). Let u G ^^(Ao,^). Then we have 

(3.15) 



!P(2)m = gQ^i{x,u)u + gl,]_{x,u)Uxl+gi^2{x,u)u^l+g2,\{x,u)u2^ 
+ {g2,2{^,u) +g2,3{x,u)) M-,1-,2 +g2,A{x,u)u2^2 



and after expansion 

2'(2)m = [^1,0,1 (■x,OOo,i+^i,i,i(.^,OOi,i+^i,i,2(.^,OOi,2] 

X [^2.0,1 (-X, OOo,l + ^2.1,1 {x, 0Ol,l + ^2,1,2 (■^, 0Ol,2] " 

= [^1,0,1^2,0,1 + ^l,l,lOi,i (^2,0,1 ) + ^1,1,201,2(^2,0,1 )] " 

+ [^1.0,1^2,1, 1+^1,1, 1^2,0.1 +^l,l,lOi, 1(^^2,1,1) +^1,1.201,2(^2.1,1) 

+ ^l,l,lOl,3(Z72,0,l)"] ";ti + [^1,0,1^2,1,2 +^1,1,2^2,0,1 +^1,1, lOi,i(Z?2,l,2) 

+ ^1.1,201,2(^2.1,2) +^1,1,201,3(^2,0,1)"] "a:2+^1.1,iOi,3(^2,1,i)"^i 

+ [^1,1,101,3(^2,1.2) +^l,l,20l,3(Z72,l,l)] Wxl"v2+^l,l,20l,3(^2,l,2)"^2 

+ ^1,1,1^2.1,1 «2xl + [^1,1,2^2,1,1+^1,1,1^2,1,2] i<;tiA-2+ ^1,1,2^2,1,2 "2x2- (3-16) 



Identifying (13.151) with (13.161 ) yields 
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Proposition 3.2. A necessary and sufficient condition for the differential operator V (2) 
defined by f lJ.iJD be decomposed into the form dJ.i^D is: 

g2,\ = ^1,1,1^2,1,1, (3.17) 

§2,2+^2,3 = ^1,1,2^2,1,1+^1,1,1^2,1,2, (3.18) 

g2,4 = ^1,1,2^2,1,2, (3.19) 

gi,i = ^1,0,1^2,1,1 +^1,1,1^2,0,1 + -t(^2,i,i) +^1,1,1^1,3(^2,0,1)", (3.20) 

gl,2 = ^1,0,1^2,1,2 +^1,1,2^2,0,1 +-^(^2,1,2) +^1,1,2^1,3(^2,0,1)", (3.21) 

go,i = ^1,0,1^2,0,1 +-^(^2,0,2), (3.22) 

= ^1,1,2£>1,3(^2,1,2), (3.23) 

= ^i,i,i£>i,3(Z.2,i,i), (3.24) 

= ^l,l,l£>l,3(^2,l,2)+^l,l,2£>l,3(^2,l,l), (3.25) 

where L = b\i\D\^i +^1,1,2^1,2- 

3.2 Factorizations of systems of nonlinear differential equations 

3.2.1 Theoretical considerations and principles 

Let A be an open subset of M" and H, an open subset of W" . Examine now the factorization 
process for systems of s-th order, (s > 2), nonlinear differential equations with n indepen- 
dent variables x= {x^ ,■•• ,x") and m>2 dependent variables u = '^ (m V " " , "'") , " = u{x) 
whose associated matrix operator, M (s) , is of the form 

^i^) = KA^P,,)]i<p,,<,n' (3-26) 

the flip^q {sp,q) are Sp,y-th order linear differential operators 

■V.<7 Pk 
'Kp,q {Sp,q) = 52 22 fp,q,k,h{^^ •,•••,• )^k,h, (3.27) 

jt=0/i=l "^ — -^~^ 

m-entries 

where fp.qxh G C (A x H, M) , Sp^q = s-l+ hp^q, ?>p,p = 1 and 5^,^ = if p / <?. 

Let A and Aq be two open subsets of M" such that Aq C A. The matrix operator (M (s) acts 

on a vector valued function u= ' (^u^ ,■■■ ,m'") G c'(^o,^) as follows 



^(^)« = [%'?(^m)]i<p,^<,„« 



£ S{p,q {Sp.q) "^ 



l<p<m 



with 

^i'-i Pk 

^,q {Sp,q) "'^ = L L //'*M {"^y ' • • •'«'") Dm ^''- (3-28) 

k=Oh=\ 

The method of factorization consists in seeking a decomposition of the matrix M {s) under 

the following form 

I 
M {s) =Y{Hi{si) (3.29) 

(=1 
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where 

Hi{si) = [T,,pA^i.p.,)],^^,,^^^ (3.30) 

and 

*/,;>.9 Pk 
Ti,p,q {Si,p,q) =Y^Yj '^'■P,'i,k,hix, ■,■■■,■ )Dtf,, (3.31) 

k=Oh=l " -^ ' 

m-entries 

mthY,[^^Si = s, Si^p^q = Si-l + 5p^ci, ahp,q.k.h G (:(Axa,R) and a,-,p,^,A:,/, G C ^^= '''"•" (A x 
a,M),/ = 2,3,--- ,Z. 

Expanding (13.291 ) leads to the relations between the unknown functions a,-,p,^ ,t,/, of 9{j{si) 
and the known functions fp,q,k,h of fW (.?) . 

3.2.2 Necessary and sufficient conditions for the factorization of systems of second 
order nonlinear ODEs 

Let A, Ao be two open subsets of R such that Aq C A, and Q. an open subset of W" . Consider 
the matrix operator 

^(2) = [%,]i<p^^<^, (3.32) 

where 

2 Pk 

'^'P = mi fp,q,k,h {x, ■,■■■,■ )^k,h (3.33) 

k=Oh=l ^ ^'"^ 

m-entries 

= fp,p,0,l{x, •,..., ^Do,! +fp,p,l,l{x, •,... ,-)Dl,l +fp,p,2,l{x, ■,---,-)^2,l 

m-entries m-entries m-entries 

and for p / ^ 

1 Pk 

^.q = X ZL fp.q.kji {X, ■,■■■,■ )'Dk,h 

k=Oh=l " -''—^ 

m-entnes 

= fp.q,0,l{x,-,...,-)Do.l+fp^q,lj{x,-,...,-)Dij (3.34) 

m-entries m-entries 

with fp^q^k,h e C{AxD.,R),x = x^. Write 5W (2) in the form 

fW-(2) = iA£i(l)-^(l), (3.35) 

where 

^•(l) = [^.M]i<;,,,<m (3-36) 

with 

1 Pk 

^i,p,p ^ 2^ 2^ ^i,p,p,k,h {X, ■;•••)■ )'-^k,h 

k=Oh=l ' ^^ ' 

m-entnes 

= ar>,p,o,i(-^' •'•••'•)Do,i+«i>,p,i,i(-^, -j---,-)!)!,! (3.37) 

m-entries m-entries 
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and fov p^q 



'Ti,p,q — ai>,^,0,l (-^J ■)•••; ■ )Do,l , 



(3.38) 



a\,p,q,k,h G (:(Axa,R) and a2,p,cj,k,h e c\Ax D.,R).Let u = ' [u\-- ■ ,m'") G ^^(Ao,^) 
Then we have 



^(2)m = [!^,,]i<p,,<„ 






(3.39) 



l<p<ni 



where 



and for py^q 



%p,p uP = fp,pfl,i (x, u) uP + fp.p^.l {x, u) < + fp.p2,\ {x, u) wf^ 



On the other hand, after expansion of ( I3.35I ). we have 



M{2)u 



'Kp,q 



l<p.q<m 



£%^ 



(3.40) 



\<p<m 



where 



'K4),pU'' — ai,p,p,l,ia2,p,p,l,l "2a: + '^1.P.P.1J X/ I-'l,ft+l('^2,/7,p,l,l)"jcMx 

h=\ 
m _ 

+ « 1 ,p,p, 1 , 1 2- 1^ 1 ,/;+ 1 ('^2,p,p,0, 1 ) Mjc '^^ + [^ 1 ,p,p,0, 1 '32,p,p, 1 , 1 + « 1 ,p,p, 1 , 1 fl^2,p,p,0, 1 






+ (^l,/),/),!,!!)!,! (fl'2,p,p,l,l ) + fl:i,p,p,l,lDi,p+i (a2,p,p,0,l ) "''] "x 



+ 



2^ ^i,p,/,o,i'32,/,p,o,i + '2i,p,p,i.iDij (a2,p,p,o.i ) 



and ior p y^ q 

%4),qU'^ = [ai,p,p,l,ia2,p,qfl,l + ai,p,q,0,l'^2,q,q,l,l + Ol,/?,/), 1,1 D 1,^+1 («2,p,?,0,l ) m"^] «? 
m _ 

+ '^l,p,p,l,l 2-Dj ^^j ('32,p,ty.0,l ) "x W^ 






+ 



2^ «l./7./,0,l«2,/,<?,0,l + ai,p,p,l.lDl,l {<^2,p,q,Q,l ) 



/=1 



Identifying ( |339l) with ( |340l ) yields 

Proposition 3.3. A necessary and sufficient condition for the differential operator fW (2) 
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defined by di. J2D be decomposed into the form ^3.35^ is: 

m 

fp,P,o,i = 2- «i,p,/,o,ifl^2,;,p,o,i +'^i,p,p,i, 1-^1,1 ('^2,p,p,o,i)) (3-41) 

1=1 

fp,p,l,l = «l,p,p,0,l«2,p,p,l,l +fl^l,p,p,l,lfl^2,p,p,0,l 

+ «i,p,p,i, 1^1,1 («2,p,p,i,i ) + <^i,p,p,i,iDi^p+i (a2,p,p,o,i ) "^5 (3-42) 

fp.p.2,i = <^i,p,p,i,i<^2,p,p,i,i, (3.43) 

= D^j^^^{a2,p,p,o,i), he{l,2,---,m}\{p}, (3-44) 

= -Dij;+i(a2,p,p,i,i), h = l,2,---,m (3.45) 

and for p ^ q 



m 



fp,q,0,l — l^<^l,p,lfi,i(^2,l,qfi,\+(^l,p,p,i.lDii{a2,p,qfi,\), (3.46) 

1=1 

fp,q.\,l = Cl\,p,p,l,lCl2,p,qfl,\ + ai,p,qX),\Cl2,q,q,\,\ + ai,p,p,l.l^l.<?+l («2,p,?,0,l ) "'^)(3-47) 

= Djj;^j(a2,p,^,o,i), /iG{l,2,---,m}\{^}. (3.48) 

3.2.3 Necessary and sufficient conditions for the factorization of systems of second 
order nonlinear PDEs with two independent variables 

Let A, Ao be two open subsets of R^ such that Aq C A, and Q. an open subset of M'". 
Consider the matrix operator 

^(2) = [%,]i<,,,<^, (3.49) 

where 

2 Pk 

^,P = L L fp^q,k,h l-^, •,•■•, • )^hh (3.50) 

^=0/1=1 "^ "^ 

m-entnes 

= fp,pfl,\{x, •,..., •)Doj +fp,p,i,\{x, •,..., •)Dij +/p,p,i,2(-^, •,---,-)Di.2 

m-entries m-entries m-entries 

+ fp,p,2,\{x, •,... ,-)D2,l +fp,p,2,2{x, •,... ,-)D2,2 

m-entries m-entiies 

+ fp,p,2,^ix, •,---,-)D2,3+/p,p,2,4(-X, •,---,-)D2,4 

m-entries Hi-entries 

and for p^ q 

1 Pk 
^.9 = jmfp,'i'k.h(.^r,---,-)'Dk,h (3.51) 

k=Qh=l ^ "^ 

m-entries 

= fp,qfl,\{x, -j-.-jODo,! + fp,q,l,l{X-, ■,■■■ ,-)Dl,l + fp,q,l,2{x, •,..., •)Di^2 

m-entries m-entries m-entries 

with fp,q,k,h eciAxQ.,M.),x= {x\x^) . Write 5Vf (2) in the form 

:^(2)=iA£i(l)-5V2(l), (3.52) 
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where 
with 



Hi{l)=[Ti,p,q],<p^q^„ 



(3.53) 



1 Pk 



Y^ 52 '^i,p,p,Kh (-^ •>•••> • )D*:,ft (3.54) 

m-entries 

<^i,p,p,Q,i{x •)••• )')Do,i +ai,p,p,\,\{x. ■,..., ■)T)i\ +aipp^i2{x •,... ,-)Di^2 



and for p^ q 



'^i,p,q —'^i,p,q,0,l{x •,...,• )Do,l , 



(3.55) 



ai,p.q,k,h G (:(Axa,M) and a2,p^cj,k,h G c\Ax D.,R).Let u = ' {u\--- ,u'") £ C^{Aq,D.). 
Then we have 



5W(2)m=[%,]i<^,^^<, 



£%<? 



q=l 



(3.56) 



l<p<m 



where 



'Kp,puP = fp,pfl,iuP + fp^p^i^iu'^^i+fp^p^i2u'^2 

+ fp,p,lA "2jf' "*" ifp,P,2.2+fp,p,2.3) U^l^l + //3,/3,2,4 "2^2 

and foxp^q 

On the other hand, after expansion of (13.521) . we have 



(3.57) 



M {2)u 



■^A 



l<p,q<m 



L^,?"'^ 
.?=i 



(3.58) 



\<p<m 



where 



'Kp^pU'' = <3L/),p,l,ia2,p,p,l,l"2x' + (^l,P,/',l,2'22,/),p,l,l +fl:Lp,p,lJ<32,p,p,1.2) "^1^2 

+ fl'l,p,p,l,2«2,p,p,l,2"2A-2 [^l,P,P,0,lfl'2,p,p.l,l +fl'l,p,p,l,ia2,p,p,0,l 

+ -^p(«2,p,p,l,l) +«l,p,p,l,lDl,p+2(fl'2.p.p,0,l)"^]";t' "*" {'^l,P,Pfi,l'^2,p,p,l,2 

+ <^\,p,p,\,2a2,p,pfi,\ +^p{<^2,p,p,l,2) + '^l,p,p,l,2^1.p+2{'^2,p,p,0,l) u'']u^2 



+ ai,p,p,\,l 2l-^l,ft+2(^2,p,p,l,l)W;t'l "^1 +fl:i,p,p,l,l 2^ ^l,/i+2('^2,p,p,l,2)W;t'i "^2 

+ fl^I,p,p,l,2 2^ Dj-^^2('^2.p,p,l,l)";c'2 W^i +fl:i,p,p,l,2 2L ^l,A+2('^2,p,p,l,2)";t'2 "v2 

m ^ m _ 

+ fl^l,p,p,l,l 22^1,ft+2('^2,p,pAl ) "x' "'^ + '^1,P,P,1,2 2^Dj jj^2('^2,p,p,0,l ) ";c'2 W^ 



A=l 






+ 



2^ '^l,p,/,0,1^2,/,p,0,l + -^p(fl'2,p,p,0,: 



/=1 
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and fov p^q 



^.qU'^ 



+ [a\,p,p,\,2<^2,p,qfi,l +a\,p,qfl,\a2,q,q,\,2 + a\,p,p,\,2^\,q+2{a2,p.qfi,l)lfl]u^2 

m _ m _ 

+ '^ 1 ,p,p, 1 , 1 2^ D J ^^2 i<^2,p,q,0. 1 ) M^i M^ + a 1 ,p,p, 1 ,2 X, ^ 1 ,ft+2 (^2,p,^,0, 1 ) "^2 "'* 









+ 



2_i '^l,p,/,0,l«2,/,?,0,l + J^p{(^2,p,q,Q,\ 



/=1 



where Xp = ai,p,p,i,iDij +ai,p,p,i,2Di,2- Identifying (13.561 ) with (13.581) yields 

Proposition 3.4. A necessary and sufficient condition for the differential operator fW" (2) 
defined by AS. 491) be decomposed into the form A3. 52^ is: 



fp,pfl.\ 
fp,p,i-,i 



fi 



p,pA,2 



fp,p,2,l 

fp,p,2,2 + fp,p,2,3 

fp,p,2A 







~ }_i'^l,p,L0,l<^2,l,p,0,l +-^p(«2,p,p,0j)i 
1=1 

= ai,p,p,0,ia2,/7,/7,l,l +ai.p.p.l,lfl'2,p,p,0,l 

+ ^p{<^2,p,p,l,\ ) + '3l,p,p,l.l^l,p+2(fl'2,p,p,0,l ) "^) 

= ai,p,p,0,l«2,/7,p,l,2 +a\,p,p,l.2Cl2,p,pfl,\ 

+ ^p («2,p,p, 1,2 ) +«l,p,p, 1,2^1, />+2(fl'2,p,p,0,l)"'', 

= ai,p,p,\,\a2,p,p,\,i, 

= «l,p,p,l,2«2,p,p,l,l +«l,p,p,l,lfl'2,p,p,l,2, 

= ai,p,p,l,2a2,p,p,l,2, 



^l,A+2(^2,p,p,0,l; 
^l,A+2('^2,p,p,l,l; 
^l,/;+2('^2,p,p,l,2; 



h€{\,2,---M\{p}, 
h = 1,2,- •• ,m, 



h 



1,2,- •• ,m 



(3.59) 

(3.60) 

(3.61) 
(3.62) 
(3.63) 
(3.64) 
(3.65) 
(3.66) 
(3.67) 



and for p ^ q 



fp,qX),l 

Jp,q,l,l 
fp,q,l,2 





2^ <^l.p.lfi,l<^2,l,q,Q,l + J^p{a2,p,q,Q,l), 



1=1 



(3.68) 



'^i,p,p,i,ifl'2,p,^,o,i +a:i,p,^,0,i'^2,,?,<?,i,i +a:i,p,p,i,i^i,£;+2(fl'2,p,^,o,i)M'^,(3-69) 

ai,p,p,l,2<^2,p,qfl,l +'^l,p,qfl,l<^2,q,q,l,2 + «l,p,p,l,2^1,?+2(fl'2,p,^,0,l ) ""^X^-VO) 



Oift+2(«2,p,?,0,i), h£{\,2,--- ,m}\{q}. 



(3.71) 
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